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1 Introduction 



Recent work has revealed unexpected connections between the functional relations and 
Bethe ansatz (BA) equations arising in integrable quantum field theory, and the be- 
haviour of certain ordinary differential equations (ODEs) In this paper we extend 
these links to cover a whole class of BA systems associated with the SU{n) Lie alge- 
bras. The papers just cited, and earlier work such as |8-11|, should be consulted for the 
background to these developments, but we begin with a brief review of the main points. 

The first example to be found concerned SU{2) BA systems at a specific value 
of the twist parameter, and associated them with certain Schrodinger equations 
The generalisation to arbitrary twist turned out to require the addition of an angular- 
momentum term to the ODE |^. In studies of the Bethe ansatz, a functional equation 
found by Baxter, the T-Q relation, is known to encode the BA equations in a particularly 
neat way; the place of this relation in the differential equations side of the story was 
found in |^. After this, the question of generalisations to SU{n) with n > 2 was 
addressed in [^,^. In |^, a reasonably complete mapping of the S'f7(3) case onto third- 
order ODEs was found, and it was suggested that higher-order equations would be the 
appropriate setting to search for the more general case. Independent work by Suzuki 
provided some more concrete support for this idea. However, this latter paper did not 
treat the completely general case (both the deformation parameter and the twists were 
restricted) and the role of a particularly important generalisation of the T-Q relation, 
sometimes called the 'dressed vacuum form', was not elucidated. For this reason, we 
return to the topic here, and show how the approach of can be extended to general n. 
One of the main new features is a way of incorporating angular-momentum type terms 
into the higher-order ODEs which encodes the BA twists in a particularly neat way. This 
simplifies the discussion considerably, and even streamlines aspects of the SU{3) case 
treated in |^. We also derive nonlinear integral equations (NLIEs) for the associated 
spectral problems, and discuss duality properties. 



2 The ordinary differential equation 

The initial ODE that one might consider has the form 

((-1)"+'^ + P{x, E))^p{x) = , E) = x"*^ - E (2.1) 

with M a positive real number, related to the deformation parameter. For n=3 this is 
the first of the equations studied in |5|, while the general-n case, but with nM restricted 
to integer values, was the subject of [Q. For each n, ( pT] ) has only one free parameter, 
namely M, and so it is clear that this equation cannot hope to incorporate SU{n) BA 
systems at general values of the twists. The results of |^,^ would tend to suggest that 
these twists should be associated with the addition of terms homogeneous with of 
the form AkX~^ 'F^!^^ , k = 2 . . . n. (It will be convenient to assume that the term with 
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k = 1 has been ehminated by a suitable redefinition of ^.) This gives n—1 further free 
parameters, exactly matching the number of twists in an SU{n) BA system. However, 
the variables {Ai^} turn out to be rather inconvenient, and it is better to take a slightly 
more indirect route. 

First, define a general homogeneous differential operator of degree one by setting 

D{s)=(^-i). (2^2) 



dx X ^ 

Useful properties of this operator are 

D{g)^ = -D{-g) , (2.3) 

D{g2 - l)D(5i) = D{g^ - l)D{g2) . (2.4) 

The first of these relates the operator to its adjoint, while the second expresses a form 
of commutativity. Now, given a vector g = {g^^gi, . . . gn-i), set 

Dig) = D{gn-i - (n-1)) L»((7„_2 - (n-2)) . . . D{gi - 1) D{go) . (2.5) 

This is a homogeneous differential operator of order n; by the 'commutativity' property 
it depends on the components {go . . . gn-i} of g in a symmetrical manner. From now 
on we also impose 

n-1 / 

En n—1) 
9^ = (2.6) 

to ensure the vanishing of the term in -D(g) proportional to x~^ dx"-^ ' Finally, we 
record one more property of -D(g) which will be useful later: its indicial polynomial (see 
for example |jl2|) is 

n-1 



/(A) = n(A-5.) . (2.7) 

i=0 

With these ingredients in place, the n^^ order ODE that we propose to study is obtained 
simply by replacing the operator in (2T) by -D(g) : 

(^{-l)^+^D{g)+P{x,E)^7P{x)=0 , P{x,E) = x""^ - E . (2.8) 



3 The fundamental system of solutions 

As in we take our cue from the approach of § to second-order ODEs. We claim 
that ( |2.8| ) has a solution y{x,E,g) such that: 

• y{x,E,g) is an entire function of {x,E,g), modulo a possible branch point at the 
origin of the complex x plane. 
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• as |x| ^ cxD in the sector | argx| < (n+l)7r/n(M+l) , 

l|~(-if^5rTi7T7^»p(-.«-V(M+l)) (3.1) 



for p = 0, 1, ... . This holds if M > l/(n— 1); more generally the WKB-like formula 
y ~ P{x, exp(— P{t, E)^/^dt) can be used in the way explained, for n=2, 

in the appendix of [Q]. The normalisation of y has been chosen for later convenience. 
• y is uniquely characterised by the asymptotic (|3.1| ) for p = and x real. 



For n=2, hMgZ, and D(g)=^ , these properties were proved in We have not 
attempted a proof for the more general case; however, see for a related discussion. 

To construct further solutions, we setj^ 

yk{x, E, g) = c^("-i)'=/2 y{io-''x, uj-'^^^E, g) , uj = e2^Vn(M+i) ^ ^3^2) 
For A: G Z, it is easily checked that yk solves (2^). If sectors Sk are defined as 

then as X — > 00 in U S^-^+i U . . . U 3^+^ ■, ( P-H ) implies that 

„(l-n+2p)M/2 

^^(_l)P^-Mi-n+2p)(A/+i)/2^^ exp(-a;-^(^^+^)x^-^+V(M+l)) . (3.4) 



We call subdominant a solution which tends to zero fastest in a given sector; up to a 
constant multiple, y^ is the unique solution to ( |2.8| ) which is subdominant in 5^. 

Wronskians between these solutions will be important in the following, so we define 

wi';:l„„k^ = W^"^^[yk„yk„...ykJ, (3.5) 

where 



(3.6) 



and fl^\x) = ^pfi{x). We will give a special status to those Wronskians whose argu- 
ments are successive integers, and in preparation for this we set 

^^k — k,k+l,k+2,...k+m-l ■ y"^-'! 





■ h /2 . 


fm 


^^^"^^[/l,/2,.../„^] =Det 


/l /2 


f 

J m 




»[m— 1] Am—l] 
L Ji J 2 


f[m-l\ 
Jm -1 



'Note that since 0;"'*^+^' = 1, the shift in E in the definition of could equally Iiave been written 
as uj"'' , so long as k remains an integer. This would have been more in line with |I|,^pT[, and emphasises 
the fact that in the 'semiclassical' limit M — > 00, the shifts tend to 1. Here we retain the conventions 
of 1^, since this makes the transition to BA equations in standard form a little simpler, and also appears 
to be the most natural choice when half integer values of k come to be discussed, in §7 below. 
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Note that = yk, and that (generahsing (|3.2|) ) 

..^-.(-'^'g) =-"^^"-"^'^'<l... kJ^-'x,^-'''E,g) . (3.8) 

For an n*^ order ODE with vanishing (n— l)*'^ order term, as is true of our case 
by (^), it is standard that all n-fold Wronskians •••/«] are independent of x, 

and vanish if and only if the solutions fi, ■ ■ ■ fn are linearly dependent. In particular, 
(|3.4D can be used as |x| — > oo in 8^^1/2 U 5^+1/2 to show that wj:"''^ = 1, and so, for 
each k G the functions {yk,yk+i^ ■ ■ ■ 1 Uk+n-i} are linearly independent, and furnish 
a basis of solutions to (|2.8D . 



4 Stokes multipliers and functional relations 

We next establish a basic functional relation. Expanding yo in the basis provided by 
{2/1, 1/2, . . . , yn}) it must be possible to write 



Y^(-l)f^C^f')(E,g)ykix,E)=0 



(4.1) 



fc=0 



where C^^\E,g) = 1 and the Stokes multipliers C^'^)(S,g), A; > are analytic functions 
of E and g. (The formal similarity between this equation and equation (5.50) of |13] 
should be noted.) Given that W^'^^ = 1, we have the simple relation 



c^''\e, 



(4.2) 



with the hat (k) indicating that the corresponding index is to be omitted. Clearly, 
C7(") = W, 



(n) 

J = 1. To treat the other Stokes multipliers, we first relate them to the 
'privileged' Wronskians W^"^\ We will only need C^^^ here, but the discussion can be 
generalised to the other cases. For general p, consider the determinants 





r m 
Vo 


yy 


m 

yp 




yo 


yi 


yp 


= Det 


y'o 


y[ .. 


y'p 






b'-i] 

yi ■■ 


[p-i] 
• yp 



(4.3) 



with i = 0, 1, . . . p—2. Expanding the first row of each using Cramer's rule, we have 

yk 







k=0 



y'k 



iVk 



w. 



(p) 



01...A;...p-l 



(4.4) 
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Now form a (p— l)x(p— 1) matrix with the RHS of ( [4. 4] ) as the first column, and the 
vectors Vj = [yj,y'j, . . .y^-' ^^]*, j = 2, 3 . . . p—\ as the remainder. Taking its determinant 
yields the following Plucker-type relation: 

Rearranging, and supplementing the notation ( |3.7D with the convention w'j^^ = 1 VA; , 



^^02...p ^ ^''02...p-l _^ 1^2 ^0 ^ '^2 '^0 u 6) 



the second equality being obtained by repeated substitution. This is the result we need. 
Since C^^^ = Wq2^ „ , we can set p = n and multiply (4^) through by 11^^=0 ^i'^ to find 

j=0 m=0 \i=0 / \j=m+2 / 



5 Bethe ansatz equations 

As it stands, equation (^]^) involves x as well as E and g, and so is not quite ready to be 
mapped onto a standard set of BA equations. The simplest idea, simply to set x to zero, 
only works at exceptional values of the gi , since in general the differential operator -D(g) 
has a (regular) singularity at a; = 0. Instead, motivated by the treatments of |2|, ^, ^, 
we expand y{x, E, g) as 

n-l 

yix, E, g) = W^'^ (x, E,g) = Y, Dil] (E, g) Xi (x, E, g) (5.1) 

1=0 



where the Xi form an alternative basis of solutions to (2.1), fixed by the demand that 



they have the simplest possible behaviours near the origin: 

Xi ~ + 0(xf'+") , x^O. (5.2) 



(Recall that the gi are the roots of the indicial polynomial ( |2.7| ).) Strictly speaking the 
asymptotic (|5.2| ) does not always suffice to pin down all of the Xi- Assume, until further 
notice, that the gj's are ordered as 

3f?e(<7o) < M9i) < ■ ■ ■< M9n-i) ■ (5.3) 

Then xo is certainly uniquely determined by (p.2|), and the Xi>o t)e defined by a pro- 
cess of analytic continuation from this solution, just as is done for the radially-symmetric 



5 



Schrodinger equation when passing between the so-called regular and irregular solutions. 
This will be used in §8 below. 

The direct substitution of into ( |3.5| ) yields an expansion for Wq"^\x, E,g) in 
which the functions Xi{^~^Xi uj~'"'^'^^E, g) appear as well as the Xii^: -^i §)• However, all 
of these functions are solutions to the initial ODE, and by considering their behaviour 
near x = one finds 

Xi(a;-^x,u;-"*"i?,g) = u;-'"^^x^ix, E,g) . (5.4) 
Therefore the m-fold Wronskians W^j^^ have expansions of the form 

W (X, ^, g) = ^InLjr.] g) ^^"^ , Xn ■ ■ ■ X^J E, g) . 

0<jl<i2<---<jm.<n-l 

(5.5) 

Using (13), (|5l|) and (U), the coefficients .^^{E, g) can be expressed as sums of 

products of the D^^) {uj~'^''^^^ E , g). This leads to relations which generalise the 'quantum 
Wronskians' of However, for current purposes it is better to treat the coefficients 
with different values of m as independent functions, and so we will leave further discus- 
sion of this point to another occasion. 

We will initially focus on the dominant terms of the expansions (^.5|). With the 
ordering ( |5.3| ), these are 

^o^™)(x,ii;,g) ~i?fo?..,n-i](^'g)^^'"^'"^"""^^'' (5-6) 
where in order to simplify subsequent calculations we set 

m—l 

Pm=Yl 9j - m(n-l)/2 . (5.7) 

Substituting (|5.5D into (|4.7| ), an x-independent equation is found by extracting the 

coefficient of the leading power x", a = X]j=o('^i • Shifting E to uj"'^'^ E 
and setting 

Z?(-) {E,g) = Z?[ot^„„_i] {E, g) , D^^^ {E, g) = {u;-^'^"'E, g) , 

T('\E,g) = C7«(u;"^-^i?,g) , (5.8) 
the final result can be written as 

n n—l /m—l \ / ™ \ 

j=0 m=0 \j=0 I \j=m+2 / 

This is a generalised T-Q relation, with the D^^^ playing the role of the eigenvalues of 
the Q (or A in ||Tl|) operators. In the language of integrable lattice models relations 
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of this type are known as the dressed vacuum form for the transfer matrix eigenvalues 
(see for example [14, 15|). 

We can now derive an initial set of BA equations. Suppose that the zeroes of 
D^"^\E) are at E = f'^\ k = 1, 2 . . . oo. (For the next few equations we will leave the 

dependence of all functions on g implicit.) Setting E = in (p.9|), the LHS vanishes, 
as do all but two terms in the sum on the RHS. Assuming, as will generically be the 
case, that there are no further vanishings, this gives us the following set of coupled 
equations for the {F^^*} : 



(5.10) 

These form a system of Bethe Ansatz equations (BAE) of SU (n) type (see, for example, 
equation (3.32) of [|l3|). However the reality properties are not very transparent when 
the equations are given in this form, and it is useful to make one more definition, setting 

Z)^'") {E, g) = A^'") (a;-"^-^(™-i)/2i?, g) , F^™) = ^-Af(™-i)/2^M _ ^^^^^^ 
Then the BAE ( |5l0| ) become 

^(m-l)(^nM/2_gM) ^(m) (^-nM£;M ) ^(m+1) (^nAf/2£;M ) 

(5.12) 

These can also be written using the Cartan matrix Cmt of the SU{n) Dynkin diagram: 

Il-''-'' \J M ^ = 1,2,.... (5.13) 

Finally, the BAE can be given a factorised form once the growth properties of the 
functions involved have been established. A WKB treatment along the lines of [Q] 
shows that the function A'^^^Ejg) = D^{E,g) has the large \E\ asymptotic 

In A^^\e, g) K{nM,n)i- FY, \E\ ^ oo , \ai:gi-E)\ < n , (5.14) 

where 

nM ' r(l + i + l)sm(f + f) ' 

The asymptotics for the remaining In A^"^^ (E, g) are more tricky. So long as it is assumed 
that the other functions d'^^^ share the asymptotic ( 5.14 ), then the quantum Wronskian- 

like relations mentioned just after ( |5.5D can be used to show that all of the A^"^\E,g) 
are of the same order, namely fj,. For this section this suffices, but for the NLIE a more 
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precise result will be required. Unfortunately, cancellations of leading terms seem to 
be at work inside the quantum Wronskians in almost all directions in the complex E 
plane, making a direct calculation rather difficult. However, indirect evidence suggests 
the following: 

ln^('")(^,g) ~ sm(7rm/n) ^ ^ \arg{-E)\ < tt . (5.16) 

sin(7r/n) 

Note that this claim matches the Z2 symmetry property of the spectral determinants 
discussed in §7 below, and its implications are in agreement with the numerical checks 
for the "soluble" M = 1/n cases reported in §8. 

For M > l/(n— 1) the order /i of the functions A^^\E,g) is less than one, so the 
Hadamard factorization theorem implies 

00 p 

^(-)(i?,g)=A(-)(0,g)n(l-n^) • (5.17) 



This finally allows ( 5.13| ) to be written as 

n-l 00 / ^(t) _ ^^c^, ^(rn) ^ 



E 




-1, A; = 1,2, ... . (5.18) 



6 The nonlinear integral equation 

Next we derive a set of coupled nonlinear integral equations. Define 



^1 ^W(cj— '^'"*S,g 



so that o(™)(£;^"'\g) = -1 by (|5l^ ). We will follow the ideas of [0, ID (another 



approach can be found in |T^). The product representation ( p.l7| ) allows Ina^™) to be 
written as an infinite sum over the zeroes now make two important 

assumptions: that all of the e'^^^ lie on the positive real axis of the complex E plane, and 
that these are the only points in some narrow strip about this axis at which a^"^\E,g) 
is equal to —1. We expect that these will hold for some range of the parameters g. 
Cauchy's theorem can then be used to replace the sum by an integral along a contour 
C which runs from +cxd to above the real axis, encircles the origin and returns to 00 
below the real axis: 

lna(-)(S,g) = — ^5]C7^iA+^ / ^^Fmt{E/E')dE>Hl+a^'\E',^)) (6.2) 



8 



where 

Fmt{E) = \n . (6.3) 

1 - a; 2 

After a variable change E = exp(^//x), we define (with a mild abuse of notation) 
lna*^™)(6') = lno('")(e^/^,g) , integrate by parts and use the property [a^'^\6)]* = 
aM(r)-i to find 

n— 1 „^ _ . n— 1 



" ^ poo O-rri 

lna("^)(e)-^ / de'Rmt(.e-e')lna^'\d') = -—— C^nt /3t (6.4) 

n— 1 „Q(3 

- 2i ^ / de'Rmt{0-e')Qm ln(l + (^'-iO)) . 

Here Rmt{G) = ■^d0Fmt{e^/^). We now take the Fourier transform of this equation, 
setting 

/oo 
fie)e-^"^d9, (6.5) 
-oo 

~ 1 /"°° ~ 

The transformed equations can be written in a compact form as 

t=l ' t=l 

-2i^^„t(A;)9mJ^[ln(l + aW)](fc) , (6.7) 
t=\ 

where the non- vanishing entries of Rmtik) are 

sinh(^^A;) ~ sinh(^(l - C)A;) _ 1 

^<'"*>^^^ = ^hn^(f(fTM ' ^™-(^) = sinh(i(i+ofe) ' ^^^M> 

the notation < mt > implying that the nodes m and t are connected on the Dynkin 
diagram of SU{n). Applying the inverse matrix (1 — R{k))~^ , transforming back and 
rewriting the imaginary parts in terms of values above and below the real axis, we 
obtain a set of coupled NLIEs for the functions a^'^\9), m = 1, . . . ,n— 1 , considered 
along contours Ci and C2 which run from —00 to +00, just below and just above the 
real ^-axis: 

lna("^)(e) = iiramig) - ih^M^e'^ 

'^J^^ilc d0'^mt{e-e')Hi + a^'\e'))- de'ip^ti9-e')Hi + -^^)y (6.9) 
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The driving terms ih^Mne^ arise from zero modes, which can be traced to the poles of 
(1 — R{k))~^ at k = i, and their magnitudes 



60 = 2 sm{TT fi) K{nM , n) , M„ 



sin(7rm/n) 
sin(7r/n) 



(6.10) 



are fixed by imposing the asymptotic ( ^.16 ). The kernel and twist factors are, respec- 
tively, 



ar(g) 



-1 



2tt 



1 - (1 - R{k)) 

n-l 



-1 



mt 



(0), 



n(M+l) 



t,m=l 



V n-R{0)) Cmtfit . 



(6.11) 
(6.12) 



These can be written more explicitly with the help of the deformed Cartan matrix 



Cmt{k) 



-1 
cosh (2^) 



m = t, 
< mt > 



and its inverse 



coth(^A;) sinh(^(n-m)A;) sinh(^t k) 
sinh(7r/c) 



(6.13) 



{m>t). (6.14) 



By ( |6.8| ), the non- vanishing off-diagonal and diagonal elements of R{k) are related by 

2R<mt>{k) = {l-Rmm{k))/cosh{f), and so 



Smt - Rmt{k) = l;C.mt{k){l - Rmm{k)) 



where 



1 Rmm{k) 



2sinh(^^A:)cosh(^fc) 



sinh(^(l + e)A:) ■ 
Plugging relations ( |6.15 ) and (6.16) into ( 6.12| ) and ( |6.11| ) we obtain 



(6.15) 
(6.16) 



dk i 



2-K 



'mt 



sinh(^(l + g)fc) 
sinh(|^A:)cosh(^A;) 



I3 

n 



(6.17) 



with ^ = 1/M. The final NLIEs exactly match the massless versions of the equations 
found in |jl8|,|l^ in a completely different context. 
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7 Symmetry properties 



The SU{n) Dynkin diagram has an evident Z2 symmetry, swapping nodes m and n—m. 
This symmetry is also respected by the BAE so long as the twists are transformed at the 
same time. We would expect to see signs of this symmetry in the differential equation, 
but quite how is not immediately obvious. For example, the first node of the Dynkin 
diagram is associated with the basic solution y{x,E,g) while its partner, the (n— 1)**^ 
node, is associated with an (n— l)-fold Wronskian Wq^ ^\x, E,g). The answer to this 
puzzle is provided by the fact that taking the Wronskian of n— 1 solutions to an n^^ 
order ODE automatically provides a solution to the adjoint equation. For n=3 this 
observation can be found in |^^, and was related to the current context in j^. For 
general n, the result can be established as follows. It is easiest not to use the factorised 



form (l2.5|) of -D(g) for the differential operator, but rather to suppose that (|2.8D has 
been rewritten as 

^V' + i?2(x)^V' + i?3(x)^V' + ... + (-ir^P(x,i?)V' = 0. (7.1) 

(For the differential operators discussed in §2, Bj.{x) = {—l)'^'^^AkX^^, but neither this 
nor the explicit relation between the ^^'s and the giS will be needed for the argument.) 
Now consider n— 1 solutions to this equation, gathered together into a vector as f = 
(/ii /21 • • • fn-i)- We denote by V[ai, 02, . . . a„_i] the determinant of the matrix whose 
rows are the derivatives ft'^il, ft"^], . . . fl''"-!]^ so that the Wronskian of /i, /2 . . . /n-i is 

[/i, /2, . . . /„_i](x, E, g) = y [0, 1, . . . , n-2] . (7.2) 

Differentiating once, and using the general property that the differential of a determinant 
is equal to a sum of determinants in which each row has been differentiated, 

— y[ai, 02, . . . a„-i] = ^ V[ai, . . . Uk-i, 0^+1,0^+1, . . . a„_i] , (7.3) 

k=l 

we have 

_1 w{n~l) ^ yrg^ ^ ^ ^ ^_4^ ^_3^ ^ (7 4) 

dx 

all other terms vanishing by antisymmetry. Differentiating again, 

VF(n-i) = yiQ^ n-4, n-2, n-ll + ^[0, 1, . . . n-4, n-3, nl . (7.5) 

The last term can be rewritten using ( [7.1] ) to substitute for f["l: 

V[0, 1, . . . n-4, n-3, n] = -B2ix)V[0, 1, . . . n-4, n-3, n-2] = -B2{x)W^''-^\ (7.6) 
(Again, all other terms vanish by antisymmetry.) Thus 

(d^ + Mx)) = V[0, 1, . . . n-4, n-2, n-1] . (7.7) 
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The pattern should now be clear. Differentiating both sides of ( |7.7D will again result in 
two terms on the RHS. The second of these can be rewritten using (|7.1| ), but this time 
yields S3(x)Ty("-i) instead of -B2{x)W^''-^^ as the one non-vanishing contribution. 
Taking this over onto the LHS and continuing to differentiate, one finally finds that 
satisfies 



which is indeed (up to an overall sign) the equation adjoint to (7.1). Now we can return 



to the factorised form of D{g), take its adjoint using (p. 3D, and rewrite (7^) as: 

-l)"+iZ)(gt) + (-l)"P(x,^))vF('^-^)(x,^,g) =0 , P{x,E) = x''^ - E , (7.9) 
with 



g 



^ = ialgl, ■ ■ ■ 9i-i) , a] = n-i - Qn-i-j ■ (7-io) 



For n even this is the original ODE ( |2.8| ), modulo the swap to the 'conjugate' set of 
parameters g^. For n odd, the term P{x,E) has been replaced by —P{x,E). However, 



the conventions adopted in the definition (3.2) mean that the 'half-shifted' functions 
Uk+i/iix^E,^ (k E Z) themselves solve (2^) with P replaced by —P. Taking an 
(n— l)-fold Wronskian of these functions and repeating the above discussion, the minus 
signs cancel. In particular, this means that whether n is even or odd we have 

-ir+'D{g^)+P{x,E))wi''_-'l(^,E,g)=0. (7.11) 

This is exactly the equation solved by the yk{x, E,g^). Comparing asymptotics as 
|x| — > oo in 5o, we finally establish the formula 

wi':_-/lix,E,g) = y{x,E,g^). (7.12) 

This shows that the (n— l)-fold Wronskian is (up to a shift in its arguments) just a basic 
solution to another ODE, and thus demonstrates that the relation between the first and 
last nodes is indeed reflected in the properties of the differential equation. The spectral 
functions corresponding to the remaining nodes on the Dynkin diagram can now be 
obtained either as m-fold Wronskians of the W^^"]^y2' ^ (n— m)-fold Wronskians of 
the Uk, and the diagram symmetry should be reflected in the equality of the results of 
the two calculations. Keeping the normalisation of the y's as in ( |3.4[) , the necessary 
identity is the following: 

Vr^":r^ = H^M [VFo^"'!) , , . . . , W^Zp] . (7.13) 

While we do not have a general proof, for m = 2 this result follows from the Jacobi 
identity (cf. for example equation (2.20) of [|l^) and the fact that VF^") = 1. This case 
is also equivalent to a formula due to Frobenius pTl , more recently reviewed in [p^] . 
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8 Analytic continuation and the linear potential 



Thus far, we have restricted attention to the dominant terms in the expansions (p. 



However, quite which term is dominant depends on the ordering (5.3) of the QiS. While 
the ODE ( |2.8| ) is insensitive to this ordering (as remarked earher, it depends on the 
Qi in a symmetrical manner), the same is not true of the NLIE (|6.9|), since it sees not 
the parameters g but rather the BA twists Oi{g). By analytic continuation, it should 
be possible to move between different twists corresponding to the same ODE, thereby 
accessing the other terms in the expansions. In this section we will test this idea by 
means of a simple example. This also serves as a rather strong check on the various 
conjectures that have been made above. 

In the expansion of y given by (|5.1|), the coefficients d\^) are explicitly 



nW^F „^ ( -I y ^^"^ [Wi^\x, E, g), xo, • • • Xi-i,Xi+i, ■ ■ ■ Xn-i] , 



where 



iy(")[Xo,...Xn-i]= n (9J-9^)■ (8.2) 

j=0,i=j+l 

Thus, for example, continuation of the parameters from g = ((70, • • • ft, • • • 5j, • • • ffn-i) 
to g' = {go, . . . gj,. . . Qi, . . . gn-i) has the following effect 

i^g) (E, g') = Z)« {E, g) , {E, g') = {E, g) , (8.3) 

with the other {E, g) remaining unchanged. So long as the property ( p.6| ) is pre- 
served at all points of the continuation, the NLIE will hold throughout. At the end of 
the continuation, the ODE will have returned to its original form, but the NLIE will 
have undergone a nontrivial monodromy, since the twists will have changed. Explicitly, 
an NLIE describing D^-^{E, g) can be found by continuing the twist parameter a to 
a' via the transformation go — > gi, with the remaining being allowed to permute 
amongst themselves. In fact, a number of different NLIEs can be obtained, depending 
on the permutation chosen. 

A set of special choices for the vector g provides a simple illustration of all this. 
We use the notation g to indicate a vector with distinct components gj taking values 
in the set {0, 1, ... (n — 1)}. There are n\ different g, but for each choice the resulting 
differential operator D{g) is equal to . We will use gi to indicate go has been fixed 
to be equal to i. In this particularly simple case, the functions (^.1[) become 



Dl}/{E,g,) = l^y{0,E,g) . (8.4) 



i\ dx 

(Since the function y{x, E, g) does not depend on the ordering of the components of 
g, we omit the subscript i on the RHS of this equation.) The direct treatment of 
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previous sections only allowed us to discuss D^^) = y{0,E,g), but according to the 
continuation idea just described, the NLIEs should also encode information about the 
higher derivatives yW. 

To make a numerical check, we set M = 1/n, so that P{x, E) = x — E and the ODE 
is solvable by a complex Fourier transform, as in ||5| for n = 3. This lies outside the 
range M > l/(n— 1), but, just as in previous examples 0,^, we assume that the final 
NLIE continues to hold. If we write y{x, E,gi) = A{x—E), then the ODE becomes 

{-ir^'^^ + xA{x)=0 (8.5) 
and has the normalised solution 

The contour T follows a path in the lower half complex plane that starts at |x| = oo along 
the ray arg(x) = — f — and ends along arg(x) = — f + (for n=3 this is illustrated 
in figure 2 of p). Evaluating ( |8.6| ) numerically for n = 3, 4 and 5, we achieved good 
agreement with results obtained via the numerical treatment of (|6.9| ), not only for the 
zeroes of the functions y{0,E,g) but also, after analytic continuation from go to gi, g2 
and so on, of y'(0, E,g), . . . , (0, E, g). The accuracy obtained was typically to 12 
digits, and could doubtless be improved with longer computing runs. Occasionally, the 
NLIE (6.9) breaks down, and in particular attempts to solve it by numerical iteration 



fail to converge. Generally this happens whenever there is a breakdown in one or other 
of the two assumptions necessary for the derivation of the NLIE in its simplest form, 
mentioned at the beginning of §6 above. For the continuations between the different 
gi it seems to be the second assumption which sometimes fails, with points appearing 
on the positive real axis of the complex E plane, not in the set {E^^}, at which the 
functions a^™^ take the value —1. (For n = 2 these points are zeroes of T{E), and the 
necessary modification of the NLIE was discussed in ||2^.) In most cases we found that 
by altering the choice of gj the problem could anyway be avoided. One exceptional 
case is n=3. The two choices for gi are (1,0,2) and (1,2,0) , and for neither of these 
did iteration of the basic NLIE converge. Despite this fact, an extrapolation allowed 
a comparison with the results of the complex Fourier transform to be made, and we 
found that the zeroes of y'(0, -B, g) were indeed reproduced. Note that in this case the 
zeroes of y{0, E, g) and y"(0, E, g) can be found from an NLIE which respects the SU{3) 
diagram symmetry |5|, but that to find the zeroes of y'{0, E, g), we are forced to use an 
integral equation which breaks this symmetry. 

9 Duality 

As in Q^^], a duality transformation can be defined which relates solutions of the ODEs 
with M > with solutions of ODEs of the same form but with — 1 < M < 0. For n = 2, 
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the case of the Schrodinger equation, this relates a confining potential to a singular 
potential. The first step is to implement a (generalised) Langer [24| transformation 

z = lnx, y{x) = e("-i)^/\(z) . (9.1) 

The use of the operators D{g) and D{g) makes this task very simple, since the relation 

d 9\ ^2^z(x^,,.J..^^ .^z( d , p-1 



dx X 



-^(-)/(z(x)) = e A + ^ - /(,) (9.2) 



allows exponentials of z to be passed successively through the factors of -D(g), resulting 
in the transformed equation 

where ji = gi — {n—l)/2. The exponentials can now be exchanged via z jj^^^ +ln 
to obtain 



where -E = —(M+1)'^^^E~^^^^^^ . Applying the inverse transformation, ( |2.8| ) becomes 
((-l)"+^Z)(g) + P(x,^))y(x) = , P{x,E) = -x''^ -E , (9.5) 

with 

~ M ^ _ _ , _ gj + M(n-l)/2 , , 

M+T ' S = {9o,9i,---,9n-i} , 9i = • 



(Observe that the property ( |2.6| ) is preserved by the transformation of g.) This, the 
dual equation, is the same as the original ODE (|2.8D apart from the reversed sign in 
front of the 'potential' term x"^, and the fact that the parameters (M, g) have been 
replaced by (M,g). Thus the equation transforms very simply under duality. 



10 Conclusions 

In this paper we have found a hidden SU{n) structure inside certain n*^-order differential 
equations. This has enabled a large set of Bethe Ansatz systems, and their associated 
nonlinear integral equations, to be related directly to spectral problems. 

One aspect that we did not discuss relates to the truncations of the functional 
relations which can occur at special values of the parameters M and g. These can 
provide a link to equations arising in the thermodynamic Bethe ansatz, a subject of 
much independent interest. The earlier analysis of [Q, which should be considered 
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as complementary to our work, was much closer to this approach, and it would be 
interesting to develop it further. 

The correspondence between Bethe ansatz systems and ordinary differential equa- 
tions has various potential applications. From the point of view of integrable systems, 
results such as duality properties are made much more apparent by the mapping into a 
differential equation, a fact which has already been exploited, for n=2, in Q. In fact, at 
least for n = 2, the T and Q-functions are coming to have an increasingly important role 



to play in the general study of integrable quantum field theories with boundaries |25-30|], 
and this makes it likely that the new perspective on these functions provided by their 
reinterpretation as spectral determinants will prove of broader relevance. (This subject 
is even starting to find applications in string theory To give one example, in the 

study of boundary flows direct relations between the T's and the (^-functions of Affleck 
and Ludwig |^2| have been found in certain cases 



On the other hand, it is also rather remarkable that one can obtain information about 
spectral determinants for quantum mechanical problems, and more general ODEs, using 



conformal field theory techniques such as the truncated conformal space approach |25] 
and perturbation theory ||2,^,3^. Our treatment also entailed some rather strong 
conjectures about the reality and positivity properties of spectral problems associated 
with the ordinary differential equations under consideration. These conjectures were 
well-supported by numerical work, but proofs are still lacking. Previous examples of 
spectral problems which, despite having no obvious hermiticity, still exhibit real spectra 
include the CV symmetric quantum mechanics of . (This was briefly discussed from 
the point of view of functional relations and T-Q systems in Q.) It appears that a 
much larger set of problems sharing similar properties is being uncovered, and it would 
be very rewarding to understand this behaviour at a deeper level. 
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